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The problem of diffusion in a time-dependent (and generally inhomogeneous) ex- 
ternal field is considered on the basis of a generalized master equation with two 
times, introduced in [1,2]. We consider the case of the quasi Fokker-Planck approx- 
imation, when the probability transition function for diffusion (PTD-function) does 
not possess a long tail in coordinate space and can be expanded as a function of 
instantaneous displacements. The more complicated case of long tails in the PTD 
will be discussed separately. We also discuss diffusion on the basis of hydrodynamic 
and kinetic equations and show the validity of the phenomenological approach. A 
new type of "collision" integral is introduced for the description of diffusion in a 
system of particles, which can transfer from a moving state to the rest state (with 
some waiting time distribution). The solution of the appropriate kinetic equation 
in the external field also confirms the phenomenological approach of the generalized 
master equation. 

I. INTRODUCTION 

Models of continuous time random walks (CTRW) [3] , for objects that may jump from one 
point to another in a generally inhomogeneous medium and which may stay in these points 
for some time before the next usually stochastic jump, are important for the solution of 
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many physical, chemical and biological problems. Recently these models have been applied 
also in economics and in social sciences (see, e.g., [4-6]). Usually the stochastic motion of 
the particles leads to a second moment of the density distribution that is linear in time 
< r 2 (i) >~ t. Such type of diffusion processes play a crucial role in plasmas, including 
dusty plasma [7] , in nuclear physics [8] , in neutral systems in various phases [9] and in many 
other problems. However, in many systems the deviation from the linear time dependence 
of the mean square displacement have been experimentally observed, in particular, under 
essentially non-equilibrium conditions or for some disordered systems. The average square 
separation of a pair of particles passively moving in a turbulent flow grows, according to 
Richardson's law, with the third power of time [10]. For diffusion typical for glasses and 
related complex systems [11] the observed time dependence is slower than linear. These two 
types of anomalous diffusion obviously are characterized as superdiffusion and subdiffusion. 

The generalized master equation for the density evolution, which describes the various 
cases of normal and anomalous diffusion has been formulated in [1,2] by introduction of the 
specific kernel function (PTD) W{r, r', r, t — r) depending on two times, which connects in 
a linear way the density distributions / of the stochastic objects (or particles) for the points 
r' at moment r and r at moment t. The approach suggested in [1,2] clearly demonstrates 
the relation between the integral approach and the fractional differentiation method [12] 
and permits to extend (in comparison with the fractional differentiation method) the class 
of sub- and superdiffusion processes, which can be successfuly described. On this basis 
different examples of super diffusive and subdiffusive processes were considered in [2] for the 
various kernels W and the mean-squared displacements have been calculated. The idea of 
the generalized master equation with two times [1,2] for diffusion in coordinate space has 
been recently used in [13] for the calculation of average displacements in the case of a time- 
dependent homogeneous external field. In [13] the jumps of the particles are assumed to be 
instantaneous, all particles are practically trapped and the electric field does not act on the 
waiting probability, which is independent of the external (electric) field. In this conditions 
the characteristic time scale of the external field has to be large (in comparison with the 
other time scales of the problem) and the probability of jumps is connected locally in time 
with the external field. As the result, in the diffusion equation the external field is placed 
outside of the integral on time. 

It should be noted, however, that in general case of the problem of diffusion in a 
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time-dependent external field the force is placed under the integral over r (see the semi- 
phenomenological consideration in [14] and Eqs. f lT5|) . ffTB"]) below). 

The general phenomenological approach to this problem has been formulated in [14]. 

This paper is motivated by the necessity to describe in more detail the influence of time- 
and space dependent external fields on the continuous-time random walks. The equation 
formulated in [1,2] is appropriate for this purpose and offers the opportunity for consideration 
of CTRW for both cases: long-tail space behavior of the PTD function, as well as for the 
fast decay of PTD function in coordinate space, when the Fokker-Planck type expansion is 
applicable. For simplicity we consider in this paper only the last case. 



II. GENERALIZED MASTER EQUATION 

Let us start from the generalized master equation with two times [1,2]: 

/(r, t) = f(r, t = 0) + ^drj dr' {W(r, r', r, t - r)/(r', r) - W(r', r,r,t- r)/(r, r)} . (1) 

Equation (CQ) can be represented in an equivalent form, more similar to the structure of the 
Fokker-Planck equation, where the initial condition is absent: 



J dr J dr' {W(r, r', r, t - r)f(r', r) - W(r', r,r,t- r)/(r, r)} . (2) 
J* dr j dr' {P(r, r', r, t - r)/(r', r) - P(r', r, r, t - r)/(r, r)} , (3) 



df(r,t) d 
dt dt 

or 

df(r,t) 
dt 

where the PTD-function P(r, r', r, t — r) is given by: 

P(r, r', r,t-r) = 2W(r', r,r,t- r)5(t - r) + ^ W(r', r,r,t- r) (4) 

Apparently, different - but equivalent - forms of the master equation exist with different 
kernels, although connected analytically. The form ()3]) is more similar to the form introduced 
first in the papers [14-16], where memory effects have been considered in a very general 
form on the basis of a master equation with one time argument t — r, which describes 
the retardation (or memory) effects. It should be stressed, that in [16], in particular, the 
straightforward connection of the generalized master equation (GME) with the usual CTRW 
model has been established. In the framework of the specific multiplicative regime of the 
function P(r, r', t — r) = P(r, r')£(£ — r) the dependence of P(r,r') and ((t — r) on the 
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waiting time distribution and the jump length distribution is quite clear (see Eqs. (9), (10) 
in [15]). The same applies to the function W, which is connected with P by Eq. (T4J). Similar 
problems for the kernel, depending on one time variable, have been discussed in [17]. In 
our further consideration we will derive the memory function as a function of the waiting 
time following the same line as in the papers [14-16] and we find the additional retardation 
function, which is the retardation of the mobility under the action of an external force 
(physically similar to dispersion of conductivity after Fourier-transformation in time). A 
description of this new retardation function depends on the specific model for the mobility 
and this will be considered in a separate paper. The argument t — r describes the retardation 
(or memory) effects, which can be connected in the particular case of multiplicative PTD 
function W(r, r', r, t — r) = W(r, r', r)x(t — T) with, for example, the probability for particles 
to stay during some time at a fixed position before moving to the next point. An equation 
with retardation, with the W function depending only on one time argument t — r, has been 
suggested first in [15] and applied in [16] to the case of the multiplicative representation of the 
PTD function. In general W is not a multiplicative function in the sense mentioned above 
and, what is more important, is a function of two times t and t—r [1]. It should be mentioned 
that the closed form of the equation for the density distribution is an approximation. In 
some cases the exact solution for density distribution can be found (see e.g. [16]- [19]), when 
a closed equation for the density distribution does not exist or gives a too rough approximate 
result. Nevertheless, in many practical situations Eqs. ([T|) or (j3j) are sufficiently exact and 
permit to describe various experimental data. 

Let us consider the role of appearance of the two time arguments in the generalized 
master equation Eq. (pQ) for the case of a time-dependent external force F(r, t). To simplify 
the consideration we can investigate the case of fast decay of the kernel W(r, r', r,t — r) = 
W(u, r, T,t — r) as a function of u = r — r', when an expansion in the spirit of Fokker-Planck 
can be applied. In this case Eq. (0Q) takes the form [1,2]: 



/(r,t) = /(r,t = 0)+ / dr 



a 



d 

A a (r, r, t - r)/(r, r) + — (B afj (r, r, t - r)f(r, r)) 



,(5) 



/o 9r a 

where the functions A a (r, r, t — r) and -B aj g(r, T,t — r)/ ff (r, r) are the functionals of the PTD 
function (the indices are equal a, (3 = x s in s-dimensional coordinate space): 

A a (r, r,t — t) = J d s uu a W(u,r,r,t - r) (6) 



5 



and 



B a/3 (r,T,t - r) = ^ J d s uu a iipW(u,r,T,t - r). (7) 

Eq. can be rewritten naturally in a form similar to Eq. fl2]), but now for the Fokker-Planck 
type approximation: 



df(r,t) d [* d 



-r / dr- 



A a (r, r, t - r)/(r, r) + J- (B Q/3 (r, r, t - r)/(r, r)) 



(8) 



<9t Jq <9r a 

We suggest, that the PTD function is independent of /(r, t), therefore the problem is linear. 

III. INFLUENCE OF THE EXTERNAL FIELDS 

One of the main sources of inhomogeneity is an external field, which also provides the 
prescribed dependence of the PTD function on r. Other words we can suggest, in the 
particular case considered, that the dependence of W(u, r, r, t — r) on the arguments r, r is 
connected with a functional dependence on the external field: 

W(u, v,r,t-r) = W(u, t - r- F(r, r)). (9) 

If an external field is absent the PTD function is a function of the modulus u = u, which 
implies that A a = and B = 5 a/ 3B (t — r) with: 

B (t - t) = / ' d s u u 2 W (u, t-r). (10) 
zs j 

For relatively weak external fields the functional (Q can be linearized as: 

W(u, t - r- F(r, r)) = W (u, t-r)+W 1 (u,t- r)(u ■ F(r, r)). (11) 

The functions Wq(u, t— r) and Wi(u,t—r) are equal to W(u, t — r; F = 0) and the functional 
derivative SW(u,t — r; F(r,r))/5(u ■ F(r, r))| F=0 respectively. Then the functions A a and 
B a p take the form 

A a (r, r, t - r) = ^F Q (r, r) J d'utfW^u, t - r) = F a {r, r)L(t - r), (12) 
where L(t — r) is given by 

L(t-r) = -l d s uu 2 W x {u,t-T). (13) 
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and 

B aj3 (T,T,t-T) = 5 a pB (t-T). (14) 



The generalized diffusion equation Eq. (JBJ) takes the form 
df(r,t) d 



j t J* dr [L{t - r) V(F(r, r)/(r, r)) + B (t - r) A/(r, r)] . 



(15) 



Of 

In general this equation contains two different functions B and L depending on the argument 
t — t. For the case of a time-independent inhomogeneous one-dimensional external field 
and in the particular case of the kernel dependence on time Lit — t) ~ (t — r) 7_1 and 
Bo(t — t) ~ (t — r) 7_1 (0 < 7 < 1) we arrive at the result, obtained in [20], [21] for the 
fractional Fokker-Planck equation. This kind of time dependence for the kernel is typical 
for the subdiffusion processes. 

The time-dependent mobility for the diffusion process (in the particular case of exponen- 
tially oscillating time-dependent external field and a time-independent diffusion coefficient) 
has been introduced in [22]. 

If the functional W{u,t — r;F(r, r)) is multiplicative, namely, W{u,t — r;F(r, r)) = 
VF(u; F(r, r))x(t — r) Eq. f JT5|) can be simplified to: 

df{r,t) d 



Jtl dTXit 



dt 

Here b and D are constants, determined by the relations: 



r) [DA/(r, r) - 6V(F(r, r)/(r, r))] , (16) 



h = -\ J d s uu 2 Wx{u) (17) 

with = 5Wiu; F(r, t))/5(u ■ F(r, r))| F=0 and 

D=7z- t d s uu 2 W {u). (18) 

As is easy to see for the external field F(r, r), which change slow in time (comparing with 
other characteristic time scales of the problem, e.g., with the time scale of the retardation 
function x{t — T )) Eq. f )16p coincides for one-dimensional case with the diffusion equation in 
[13]. 

The physical meaning of the multiplicative structure of the functional W is that the inde- 
pendence of the time delay of the random walkers is independent of the external field. The 
dimensionless function x(t) i n this simple case is associated with the hopping-distribution 
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function ip(t) = \ip*(\t) introduced in the master equation by Scher and Montroll [15], 
with A = 1/Vq (t is the characteristic waiting time for the hopping-distribution). Laplace 
transformations of these functions xi z ) an d ip*{z) relate them as follows 

«« - M-y w) 

For an exponential hopping-time distribution if)(t) = Xexp(-Xt), where A = 1/tq we have 
ip*(z) = 1/(1 + z), x( z ) — I/- 2 an d x{t) = x(^) = 1- 111 this case Eq. f|T6|) reduces to the 
usual diffusion equation in an external field with diffusion coefficient D and mobility b: 

= DA/(r, t) - 6V (F(r, f)/(r, t)) • (20) 

IV. HYDRODYNAMIC APPROACH 

In order to better understand the situation on the basis of a non-phenomenological ap- 
proach, let us consider the charged particles with an inhomogeneous density in the external 
electrical field in the hydrodynamic approximation. The equation for the density n(x,t) 
reads 

d 

—n{x,t)+divj{x,t) = 0, (21) 

where j(x, t) = n(x, t)v(x, t) and v(x, t) is the hydrodynamic velocity. In the hydrodynamic 
approximation, when the charged particles particles (with charge e and mass m) move in the 
medium under the action of an external time-dependent electrical field E(x, t) the equation 
of motion has (for constant temperature T) the form 

d 

— [n(x, t)vi(x, t)] + V fc [n(x, t)vi(x, t)v k (x, t)] (22) 

T e 

= V«n(x, t) H EAx, t)n(x, t) — vn(x, t)vAx, t). 

m m 

Here v is the effective frequency of collision with the particles of the thermostat. In the 
linear by v approximation the solution of Eq. fl23|) gives the closed expression for the flux j 
via the density n(x,t). This solution for time- independent v has the form 

j(x,t)= / dt'exp\-u(t -t')] (— [n(x,t')E(x,t / )] - — Vn(x,t')\ . (23) 
J ~co [m m J 

Inserting this value of j(x, t) in Eq. (121]) leads to the diffusion equation 
dn(x, t) r ' 



dt 



I dt' {D(t-t')An(x,t') - efi(t - t')V[n{x,t')E{x,t')]} 

J — oo 



(24) 
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where in the case considered the "effective diffusion function" and "effective mobility func- 
tion" are given by D(t) = T exp(—ut)/m and fi(t) = exp{—vt)/m, respectively. If the func- 
tions E(x, t) and n(x, t) change in time very slowly (the characteristic time for its change 
r 3> 1/u) Eq. ( 1241) reduces to the standard form of the diffusion equation 

dn ^ l) = D An(x, t) - e/ioV[n(x, t)E(x, t)]. (25) 

Here we introduced the notations D = T jmv for the diffusion coefficient and /^o = 1/mu 
for the mobility coefficient. 

Equation f l24l) represents a particular case (in hydrodynamic approximation) of the gen- 
eral relations between the fluxes and acting thermodynamical and the external forces. Of 
cause, the time integration in Eq. f[2~4"l) can be considered in the normal hydrodynamical 
conditions as an excess of accuracy due to the inequality r ^> 1/u. For us, however, the 
most important result is the general structure of Eq. (!24j) . which demonstrates that the time 
integral includes the electrical field E(x, t). The structure of Eq. f|24|) confirms the result of 
our consideration on the basis of the generalized master equation for diffusion [14], where 
the time dependent electric field included in the time integration. 

Since the equilibrium density in the external time- independent potential (p(x) has a form 
of the Boltzmann distribution n(x) ~ exp [— (p(x)/T], the diffusion and mobility coefficients 
satisfy the Einstein relation D = fi^T. In the considered case the same statement is valid 
also for the effective diffusion and mobility functions D(t) and n(t), namely D(t) = Tfi(t). 
The general structure of the diffusion equation ( |24|) is similar to the phenomenological 
Eq. ( Ti~6l) (with the appropriate renormalization of the kernel, which eliminates the external 
derivative of the time integral). 



V. KINETIC APPROACH 



Let us start with the kinetic equation for the distribution function in an electric field 

df{p,x,t) df(p,x,t) df{p,x,t) 

dt + V TTx + eE ( x ^) Qp =Ist{p,x,t). (26) 

Here I s t is some kind of "collision integral" , which can describe in general, as we show below, 
not only real collisions of particles, but also (for the appropriate problems, e.g. moving of 
the alive objects) the more complicated processes, as the displacements with some pauses, 
etc. 
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For simplicity we consider the one-dimensional case s=l, but the generalization for the 
cases s=2,3 is trivial. The distribution function f(p, x, t) is normalized to the density n(x,t) 

J dpf(p,x,t) = n(x,t). (27) 

For the case when the collision integral conserves the total number of particles, i.e. 

dpl st (p,x,t) = 0, (28) 



integration by p leads to the continuity equation 

dn(x, t) 



dt 



+ divj(x,t) = 0. (29) 



To calculate the flux j(x,t) let us use the Fokker-Planck approximation for the collision 
integral I s t(p,x,t) and rewrite for this case Eq. ( J26l) in the form 

df( P , x, t) + ^jij^ + eE{x ^ t) djij^ = a f X) t) + m 2 D m^\ m 



dt dx dp dp V ®P 

We suggest that the friction (3 and the diffusion D coefficients in velocity space are the 
constants, which satisfies the Einstein relation (3T = mD. Integrating Eq. ( 1301) by p leads 
to the expression 

dj(x 1 t) d 



dt dx 



dpv f(p,x,t) 



- —E(x, t)n{x, t) = -pj(x, t) (31) 



If we assume that f(p,x,t) has the quasi-equilibrium form f(p,x,t) = n(x,t)f (p), then we 
arrive at the following solution of Eq. ( 13TT) similar to ( 1231) 

j(x,t)= [ dt'exp[-f3(t-t')]{ — [n{x,t')E{x,t')]-<v 2 >Vn(x,t')\, (32) 
J-oc lm J 

where for the Maxwellian distribution fo{p) = Fm{p) in one- dimensional (s = 1) case 
< v 2 >= T/m. In this case the diffusion equation is equivalent to Eq. ( 12^1) obtained in the 
hydrodynamic approach, but with the change v — > (3 in the functions D(t), n(t), as well 
as in the coefficients Dq and fio- The function D(t) is naturally connected with the time 
dependent conductivity cr(t) = e 2 n Q fi(t), where n is the average density of the particles. In 
the simple case considered the respective frequency- dependent conductivity u(oj) is 
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Let us now consider the alternative case of the kinetic equation (1261) . when the collisions 
are negligible (I st = —ef(p, x, t) with e — > 0). We also suppose that the electric field is weak 
and can be considered as a perturbation. To find the evolution of the density we split the 
distribution function in two parts: f(p, x, t) = fo(p, x, t) + fi(p, x, t), where the perturbation 
fi is proportional to the electric field E(x,t). The respective kinetic equations are 

+ v 3_mj^a = . A = Mx _ vUp) (34) 

df 1 <^ + v df 1 (i^ + (35) 

The continuity equations follow from Eqs. (135|) : 

^^- + divj o {p,x,t) = 0, (36) 

where jo(x,t) describes the flux without the electrical field, and 

^^- + divj 1 (x,t)^0, (37) 

where ji(x,t) describes the perturbation of the flux in the lowest order of the electric field. 
The solution of Eq. ( |35l) reads 

f 1 ( P ,x,t) = -e [ dt'exp [-e(t - 1')} ~ Vt,P ^ E{x - v(t - t'),t'). (38) 

Now we can calculate j(x, t) = jo(x, t) + j±(x, t): 

jo{x,t) = J dpvf (x - vt,p) (39) 

ji(x,t) = / dpvfi(p,x,t) = -e / dt'exp [—e(t — t')] x (40) 

J J — oo 

[dpv dfo{x : Vt > P) E(x-v(t-t>),t>) 
J op 

The latter equation can be rewritten as 

j 1 (x,t) = -e f dt'exp[-e(t-t')]x (41) 

J — oo 

dx > [ d pv dfo{ ~ PlX '~ Vt ' ) 5{x - x' - v(t - t'))E(x', 
J op 

= I dt'exp[-e{t -if)} dx'ir(x,x' ,t,t')E(x' ,t'). (42) 

J —oo J 
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In equation Eq. (j4"2j) the function ir(x,x',t,t') is equal to 

dfo(x -vt,p) 



Tr(x,x',t,t) — —el dpv 



dp 



5(x -x' - v(t - If)), 



(43) 



in which fo(x — vt,p) can also be written as fo(x' — vt',p). The function ir(x,x',t,t') takes 
into account the processes of space and time dispersion for the inhomogeneous and time- 
dependent distribution / = fo(x — vt,p). 

Let us choose the distribution function f in the natural form f (x — vt,p) = n (x — 
vt)fo(p). Then finally we arrive at the expressions for the fluxes jo(x,t) and ji(x,t): 



jo(x,t) = / dpvn (x-vt)f (p) 



(44) 



j^x, t ) = -e J dt'exp [s(t - t')] j dx' j dpv ^^^ ^ 5(x - x' - v(t - t'))E(x', t'j]45) 



dp 



The expression for n(x, x', t, t') can be rewritten in the form 



tc(x, x' ,t,t') — —el dpv 



n (x - vt) - Ay ( p )Vn (x - vt) 

op m 



8{x - x' - v{t - 0146) 



Here and in what follows the operator V x acts only on the function no placed behind it. 
After integration by v we find 



7l(x, x' t, t') 



-em 



x — x 
Jt~V)'- 



m 



\( U 4-1 \ 114- ./M^/0(P) I 

n [(x t-xt)/{t-t )\ — \ p =m(x-x')/(t-t')- 
Vn [(x't - Xt')/(t - tf))fo{p)\p=m(x-x')/[t-t') 



(47) 



If E(x, t) is an oscillating function proportional to sin(ut) or cos(ut) or a function damp- 
ing in time, the argument x — vt under the integral in Eq. (l4"5j) equals to (x't — xt') /(t — t'). 
The expression of the particle density n (x — vt) (due to the presence of 5-function) in the 
limit of large t can be taken equal to x'. In this case the function 7r can then in good 
approximation be written in the form 



tt(x, x' , t, t') = —em 



(t - V) 



n (x 



,^/bb) 
dp 



\p=m{x-x')/(t-t')- 



t t\ 
— fo(p)\p=m(x-x')/(t-t')^ 'xno[x' '(1 + — 

m v " v ' t 

Therefore, the current ji(x,t) for large t takes the form 



t' 



x- 



(48) 



ji(x,t) 



op m 
5(x - x' - v(t - t'))E(x', t')(49) 
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Then we arrive at the approximate expression of the "hydrodynamic" electrical flux in 
the collisionless case: 



ji(x,t) = — e J dt'exp [— e(t — t')} J dx' j dpv x 
t> 

no(x')fi'(x — x',t — t') H V x ino(x')/j,"(x — x',t — t' 



m 



E(x',t 



(50) 



where the generalized mobilities are given by 

H'(x,t) = - 



dpv ^-^— 5{x — vt) 
op 



(51) 



and 



dpvfo(p)8(x — vt). 



(52) 



We can also introduce the mobility operator p, 

ji(x,t) = e / dt' I dx'exp [— e(t — t')]E(x', t')jx{x, x' , t, t')n (x') 

J — oo J 

where jl(x, x', t, t') equals 
fj,(x, x', t, t) = - 



(53) 



dpv6(x — x — v(t — t')) 



df (p),. ( ,t' 
op m 



(54) 



Therefore, equation ( 1371) for the flux perturbation associated with the presence of the 
weak electrical field in the collisionless limit has the form 

ft 



drii(x, t) 
dt 



+ 



eV x j dt' j dx'exp[-e(t -t')]E(x',t') j dpvjl(x, x' ,t,t')n (x') = 0. (55) 



If the space dispersion is negligible jl(x, x', t, t') ~ 5{x — x 1 ) and fl55l) transforms into 

f / dt'exp [-e(t - t')]fl(t, t')V x [E(x, t')n (x)} = 0. (56) 



dn 1 (x,t) 



dt 

Finally, for the case of slow changing in space of the density profile no(x), when the pa- 
rameter To < v > /I < 1 (<o > , ro and L are the average velocity of the particles, 
the characteristic time scale for the electric field and the characteristic space scale for the 
density n (x) respectively) the second term in brackets Eq. ( 1541 can be omitted and the 
operator fx modifies to the function (15T1) /i'(x — x', t — t'): 



p,(x — x',t — t') — > n'{x — x',t — t') 



dpv8(x-x' -v(t-t'))^^. (57) 
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Then the diffusion equation ( 155]) simplifies to the form typical for the case with an electric 
field present: 

dni ^ + eV x y* tft 7 y dx'exp [s{t - t')]E(x', t')^{x -x',t- t')n {x') = 0, (58) 

Evidently the function n'(x, t) is simply connected with the conductivity a(x,t) (in the case 
considered with the collisionless conductivity) by the equality a(x,t) = e 2 rio(x)fj,'(x,t). 

This consideration provides the evident answer on how the time-dependent electrical field 
should be included in the diffusion equation and permits to make the choice between the 
different forms of the diffusion equations considered earlier [14]. The structure of Eqs. (124"]) . 
(132]) and fl58|) confirms the result of the generalized diffusion equation, introduced in the pa- 
pers [3,4] (on the example of some particular form of the kernel in the kinetic approximation 
considered above). 

VI. STOP-MOVE COLLISIONS 

Now let us consider on the kinetic level the problem of transport for the particles, which 
can move in a time-dependent external electric filed as the quasi-free particles, but can be 
trapped and stay in the rest state during some time. The similar problem has been consider 
for the time-independent external field on the basis of generalized Fokker-Planck equation 
in [23]. 

Let us introduce a "collision" integral /, that takes into account the specific "jumps" of 
the particles: 

I = -uf(p,x,t) + v [ dt'tfj(t-t')f(p,x,t'). (59) 

J t 

Therefore the kinetic equation reads 

df(p,x,t) ^ df(p,x,t) df(p,x,t) = 

dt dx 'dp 

-vf(p,x,t)+v f dt'^(t-t')f(p,x,t'). (60) 

J t 

This "stop- move" collision integral describes the moving particles, which may change 
from a "moving" state to the "rest" state and vise versa. We assume that the change from 
the "rest" state to "moving" state takes place with recovering of the momentum distribu- 
tion. The momentum distribution of the moving particles, which leave the phase volume 
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{dx, dp} at the moment t' at the point of the phase space x, p is equivalent to the momentum 

distribution of the particles, which arises from the "rest" state at the position x for t > t', 

with the delay time t — t'. More complicated situations will be considered in a separate 

study. The function ip(t) characterizes the probability for the particles to stay in a state of 

rest during a time span t — t' . 

Let us consider the conservation laws for the kinetic equation with such jumps. The 

continuity equation reads 

dn f (x,t) f tiii 

— h dn;j(x, t) = / dpl(p,x,t) = —vnf(x,t) + v / dt'ip{t — t')rif(x,t'). (61) 

^ J Jt 

We have distinguished between the "flying" particles and the particles at "rest" state. 

The function f(p, x, t) is the distribution of the "flying" particles (p ^ 0). We also introduce 

the density of the "rest" (p = 0) particles n r (x,t). We use the "stop-move collision" term 

for the process of transferring between the "flying" and the "rest" states. 

The conservation of the total number of particles reads 

dx[rif(x, t) + n r (x, t)] = N, 

N = N f + N r , (62) 

where N is the constant. There is also the evident equality 

() " r{ ''- f) im f (x,t)-v f dt'ip(t -t')n f (x,t'). (63) 



dt 

From Eqs. (jBOjl. flES) it follows that 

dn r (x,t) dnf(x,t) 



t 



dt dt 

Equations for the numbers of "free" and "rest" particles are 

dN f (t) 



+ diuj(x,t) = 0. (64) 



dt 



-vN f (t)+u [ dt'i)(t-t')N f (t'), (65) 



9Nr ® -- uN f (t) - v [ dt'^{t - t')N f {x, t'). (66) 



to 



dt 

Integration of Eq. fl6~4l) by x leads to Eq. 

Now let us integrate the kinetic equation by p with the multiplier p. The relevant equation 
of motion reads (dimension s = 1) 

+ f dpv ^^ ~ e -^n f (x,t) = 



dt J dx m 

-vj(x, t) + v [ dt'ip(t - t')j(x, t'). (67) 

Jto 
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We will assume that the integral term with f(p,x,t) in Eq. (167]) can be represented as 
d(t) dnf(x,t)/dx. This the exact representation is exact for, e.g., such a form of the distri- 
bution function f(p,x,t) = f(p,t)nf(x,t). The function d(t) in this case equals: 



d{t) 



dpv 2 f(p,t) 



(68) 



For the Maxwellian distribution d(t) is time independent d(t) = d = T/m, where T is the 
temperature. In general d(t) =< v 2 > is the average velocity of the "flying" particles. 
Then Eq. (ISTj) represents the integro-differential connection of j(x,t) and rif(x,t): 



9j(x,t) , dn f (x,t) eE(x,t) . 

d\t) — — nf[x,t) 



dt 



dx 



rn 



-vj(x, t)+u dt'ip(t - t')j(x, t'). 

J to 



(69) 



In order to solve this equation we use the adiabatic switched process for "hopping collisions" 
(to = — oo) and the Fourier-transform of Eq. (1691) by time: 



where 



and we denote 



{-iu + v[l - ip(u)}}j(x,u) = <p(x,u), 



i/j(lu) = I dTexp(iuT)i/j(T) 
'q 



, .drif(x,t) eE(x,t) , . 
ip(x,t) = -d(t) — + — ) z ±-t-n f (x,t). 



The solution for the flux is then 



dx m 



duj exp(-iut) 



or 



, f du exp(\—iu(t — t')] r /s dnf(x, t') eE(x.t') . 



2tt iu — v[l — ip{io)\ dx in 

The flux can be rewritten by introducing the function %(t — t') 

, dnf(x,f) eE(x,t') 



j(x,t) = / dt'x{t-t') 



d(t> 



dx 



rn 



Tlf(x, t') 



where 



duj exp[—iuj(t — t')} 
2ni lo + iv[l — i/)(u))] 



(70) 



(71) 



(72) 



(73) 



(74) 



(75) 



(76) 
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Inserting this flux into the continuity equation we find the diffusion equation in the form 
dnfi ^ t] = J dt' X (t-t') {d(t')An f (x,t') - ^V[E(x,t')n f (x, *')]}, (77) 

which, for time- independent d, is the particular case of Eq. (124T) . based on the general master 
equation for diffusion, introduced in [3,4]. An essential feature of the diffusion process is the 
character of the influence of the time-dependent external field placed in Eq. (IT?]) under the 
time integral. This equation coincides formally with the hydrodynamic equation Eq. ( 1241) if 
x(t — t') is the retarded function — t') = for t < t'). 



VII. CONCLUSIONS 



We show that the generalized master equation with two times, which has been introduced 
in [1,2] and [14], can describe the influence of inhomogeneous and time-dependent external 
fields on the diffusion processes. Linearization of the general master equation in the external 
field leads to essential simplifications. In this case the diffusion processes depend, in gen- 
eral, on two different functions of time, which describe retardation, or frequency-dependent 
mobility and diffusion, in particular, due to the finite time of occupation and transferring 
particles in space in the presence of the external field. Relations with simpler models are 
established. The rigorous consideration on the basis of the hydrodynamic approach and 
various kinetic equations confirms the results of the phenomenological approach of the gen- 
eralized master equation. Of cause, the kernel functions W or P can only be defined in a 
concrete way in the framework of particular physical models, e.g., on the basis of kinetic 
theory with specific collision integrals, describing the stochastic motion with retardation. 
We also introduced the new stop-move collision integral, which describes the processes of 
diffusion with particles continuously changing from moving to resting and back. The ap- 
propriate kinetic equation is solved for a time-dependent external field, which also confirms 
the results of the diffusion master equation approach. This type of motion is very common 
in Nature and the introduced collision integral can easily be generalized to more complex 
processes of " stop- move" motion. The analysis presented in this paper opens opportunities 
to consider a wide class of the problems of normal and anomalous transport in external 
fields on the basis of the generalized master equation with two times. The Einstein relations 
in general are not applicable to the case of the non-stationary external field, but in the 
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particular case can be valid for the time dependent diffusion and mobility functions, as it 
was found above in the present paper. 
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